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The dynamic conductivity (j(uo) of strongly correlated electrons in a symmetry broken phase is 
investigated in the present work. The model considered consists of spinless fermions with repulsive 
interaction on a simple cubic lattice. The investigated symmetry broken phase is the charge density 
wave (CDW) with wave vector Q — (7r,7r,7r)^ which occurs at half-filling. The calculations are 
based on the high dimensional approach, i.e. an expansion in the inverse dimension 1/d is used. 
The finite dimensionality is accounted for by the inclusion of linear terms in 1/d and the true finite 
dimensional DOS. Special care is paid to the setup of a conserving approximation in the sense of 
Baym/Kadanofi' without inconsistencies. The resulting Bethe-Salpeter equation is solved for the 
dynamic conductivity in the non symmetry broken and in the symmetry broken phase (AB-CDW). 
The dc-conductivity is reduced drastically in the CDW. Yet it does not vanish in the limit T ^ 
Oh! due to a subtle cancellation of diverging mobility and vanishing DOS. In the dynamic conductivity 

cr(aj) the energy gap induced by the symmetry breaking is clearly discernible. In addition, the vertex 
corrections of order 1/d lead to an excitonic resonance lying within the gap. 
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I. INTRODUCTION 
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^\ . The investigation of the transport properties of highly correlated fermionic systems has attracted much attention 
in recent years. A thorough understanding of the conductivity in particular is essential for the technical application 
of materials such as metallic oxides in electronic devices. The development of a new analytic approach, the limit 
of infinite dimension for fermionic systemsu'D, allowed the numerical description of the metal- insulator occuring-iii 
the half-filled Hubbard model in d = cx) for higher values of the interaction U assuming a homogeneous phaseHu. 
The latter assumption means that one deliberately ignores the possible occurence of symmetry breaking for the sake 
C ' of simplicity. It is argued that on frustrated lattices symmetry breaking is suppressed so that the metal-insulator 

= transition occurs at higher temperatures than those at which symmetry breaking sets in. 
^ ' With this background in mind, it is the aim of this work to extend and to complement the results known so far 
O ■ into two directions. First, the finite dimensionality of realistic systems, i.e. mostly d = 3, shall be included at least 
O • to lowest non-trivial order in an expansion in 1/d. Much care is used in including these correction without physical 
^ , and/or analytic inconsistencies. It is shown that it is not sufficient to use a conserving, ^-derivable approximation 
JJ"! ■ in the sense of Baym/Kadanoff. Furthermore, the true three-dimensional DOS will be used. Second, the influence 
r> ' of symmetry breaking on the conductivity, especially the question of possible metal-insulator transitions induced by 

. symmetry breaking shall be investigated. 
" " ' To this end, the model of spinless fermions with repulsive interaction for particles on adjacent sites is considered 
on a generic bipartite lattice, namely the simple cubic lattice. Its Hamiltonian at half-filling n = 1/2 reads 



^ <ij> <i,j> 



where cf (c^ ) creates (annihilates) a fermion at site i. The sum X]<i j> I'uns over all sites i and j which are nearest 
neighbors. The coordination number_Z = 2d = 6 appears for the proper scaling of the kinetic energyOI and for the 
proper scaling of the potential cnergyO. The interaction constant is U. r-. 

In this model the symmetry is broken yielding an AB-CDW at half-fillingS for infinitesimal values of the interaction 
at T = and for sufficiently large interaction at all finite temperatures. The AB-CDW consists of alternating sites 
with a particle density above (below) average. The order parameter b is the absolute deviation of the particle density 
from its averageQ. As far as the occurence of a symmetry broken phase is concerned, the model of spinless fermions 
at half-filling is similar to the Hubbard model at half-filling which displays antiferromagnetic behavior. The main 
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differences are that the broken symmetry for spinless fermions is discrete whereas it is continuous in the Hubbard 
model, and the fact that a local interaction like the one in the Hubbard model does not favor a spatial order by itself. 
The latter fact leads to a value of Tc oc 1/U for large U in the Hubbard model whereas one has Tc cx J7 in the spinless 
fermions model. 

The article is organized as follows. Succeeding this introduction it is discussed how a thermodynamically and 
analytically consistent extension of the limit Z — > oo can be performed. Next the basic equations for the extension 
to linear order 1/Z are derived and their numerical evaluation is sketched. This third section contains also results 
for the DOS and the corresponding proper self-energy. In sect. 4 the Bethe-Salpeter equation is set up and solved 
for the conductivity (t{uj). The preservation of the f-sum rule is discussed. Numerical results for the dc- and the 
ac-conductivity are presented in sect. 5. The findings are summarized and dicussed in the final section. 

All energies (temperatures, respectively) throughout this article will be given in units of the root-mean-square of 
the "free", i.e. non-interacting, density-of-states of the lattice model concerned. All conductivities will be given in 
units of e^/{ha'^~'^) where a is the lattice constant. The constants a, ?i, and fee (Boltzmann's constant) are set to 
unity. 

II. PROPER SELF-CONSISTENT EXTENSION 0¥ Z = oo 

In the case Z — oo, the evaluation of diagrams and the treatment of quantities like the DOS is conceptually simple. 
It is always the leading contribution in 1/Z and only this which must be kept. There is no dependence on the sequence 
in which certain quantities and the equations relating them are considered. All sum rules which hold in any dimension 
also hold at Z — oo, continuity provided for the limit Z ^ oo. This simplicity is lost as soon as corrections in 1/Z 
are to be included. For concreteness, let us consider the linear corrections 1/Z; the problems are illustrated for the 
free DOS, the Dyson equation and the free energy F as function of the order parameter h. 

The DOS is a non-negative function of which the zeroth moment is unity. This holds in any dimension, hence in 
Z = 00. On including the linear correctional one realizes that the approximate expression becomes negative at large 
values of a;. This is a disadvantage of the otherwise systematic expansion. Another inconvenience catches the eye in 
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FIG. 1. Non-interacting DOS in d = oo (short-dashed curve), in d = 3 (sohd curve) and the DOS expanded in 1/d evaluated 
in d = 3 (long-dashed curve). These densities of states are symmetric about the y-axis. 

The expanded DOS does not improve considerably the agreement with the true finite dimensional DOS (here d — 2>). 
A finite expansion in 1/Z cannot produce the van-Hove-singularities. 

To circumvent the problem of the DOS expansion, we decide to use the exacLfinite dimensional DOS, i.e. the d = 3 
DOS. This procedure provides often even in c? = 1 a remarkable agreementl'Ll. In d = 3, this appcpximation yields 
qualitatively agreement for the local DOS as compared to finite dimensional perturbation resultaS. Presently, the 
approach of using a finite dimensional DOS in an otherwise infinite dimensional calculation as approximation for the 
finite dimensional problem is employed as so-called "dynamical mean- field theory" □ or "local impurity self-consistent 
approximation"^. 

Next the problem of a systematic 1 /Z-expansion is discussed for the Dyson equation. It is stated in a simple case 
when the self-energy is strictly local in real space, i.e. constant in momentum space 
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gH=5o(^-SH). (2) 

This case is realized, for instance, in the Hubbard model in d = oJi. No lattice site or spin index appears since the 
phase is assumed to be homogeneous and non-magnetic. The quantity g{uj) stands for the full local Green function 
Gi^i{ijj) and 5o(^) stands for the free Green function G'o;i,i(cLi). The expansion of the Greeji function corresponds 
to the expansion of the thermodynamic potential since they depend linearly on each otherll3. An expansion of the 
self-energy, however, yields a different expression for (7(0;) since g^iiio) is not a linear function. The expansion of 
the self-energy seems more promising since it preserves the Dyson equation by construction. Moreover, it is able to 
describe the shift of singularities, e.g. the shifts of the band edges. (Note that we discuss here finite expansions of 
the quantities considered). 

In spite of the choice to expand the self-energy some ambiguity persists. 
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FIG. 2. Externally applied field /iext as function of the order parameter h for (7 = 9 and T = in d = 3. The short-dashed 
curve depicts the 1/d self-consistent result, the long-dashed curve the result of a systematic expansion of the self-energy. The 
zeros of the curves correspond to thermodynamic equilibrium. But only zeros with positive slope are locally stable (6 « 0.48). 

In fig. ^, this problem is illustrated. It arises in the description of spontaneous symmetry breaking. Two results 
for the dependence of the conjugated field on the order parameter are opposed. The data refers to the AB-CDW 
occuring in the spinless fermion problem at half-filling. The dotted curve results from a fully self-consistent calculation 
whereas the dashed curve results from a systematic expansion of the self-energy. Note that the self-consistent approach 
generates higher order contributions. 

The argument results now from the strange behavior of the dashed curve in the vicinity of the origin. The free 
energy belonging to the dotted curve can be found by integration; it has an unstable maximum {dh^xT / db < 0) at 
6 = and two stable minima (9/iext/96 > 0) at 5 « ±0.48. But there is no free energy belonging to the dashed curve 
since it would have three maxima in sequence around b ~ which is mathematically impossible (theorem of RoUe). 
This is a very strong argument in favor of a self-consistent calculation. 

For completeness, it shall be mentioned that one may argue that in the vicinity of the physical solutions, i.e. the 
minima, the difference of both approaches is negligible. There are also cases known where the systematic, non self- 
consistent approach yields better resultsEa. But there is still anothespadvantage of the self-consistent treatment which 
will be crucial for what follows. In the sense of Baym/KadanofOEJ it covers also the calculations of two-particle 
properties and ensures the preservation of sum rules. So, Schweitzer and CzychoU resorted, in their calulation of 
resistance and thermopower for the periodic Anderson model to the self-consistent treatmentllJ although their results 
for the local DOS did not necessarily favor this approachO. 

As result of the above discussiDp_^the starting point for the inclusion of 1 /Z correction is the generating functional 

according to Baym/KadanofOEj. This is the quantity which is expanded in a 1/Z series. Then the truncation of 
this series yields an jSipproximation to the corresponding order. The power counting for the diagrams of <I> has been 
explained previouslyi3ii3. 
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FIG. 5. The self-energy diagrams derived from fig. ^ by taking out one propagator line. The diagrams shown contribute in 
order IjZ. 

Note that the Fock diagram is seemingly of another order, namely 0(1/Z'^/^), than the third diagram, 0(1/^), 
which is called the local correlation diagram henceforth. What matters, however, is the order relative to the free 
Green function which is 0(1/Z^/^) for adjacent sites. It is another advantage of the Baym/Kadanoff formalism that 
one does not need to bother about these questions once the approximate ^-potential is chosen. 

Now a point shall be highlightened which has not been mentioned before to our knowledge. In spite of the many 
arguments in favor of the Baym/Kadanoff formalism its naive application does not guarantee the absence of unphysical 
results. A counter example serves as illustration. Consider an approximate consisting only of the diagram in fig. 
^(a), summed over all sites z, j, such that i and j {i' and /, respectively) are adjacent to one another and fulfill 
2 ^ j' and i ^ i' . The resulting nearest-neighbor self-energy Ej j = = t' (uj) ~\- it" [uj) has a finite imaginary part 
t"{(jj). Using the Dyson equation, one obtains in the homogeneous phase 
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_ oj-{l + t')e{k) + it"e{k) 

~ (cc; + Oi - (1 + i')e(k))' + (<"£(k))' ' 

By choosing an appropriate wave vector k at fixed w one can have ilic sign of e(k) such that the imaginary part of 
Gk,k(^ + Oi) is positivetS. This is a contradiction to the exact resultt2l. Note that the details of t{uS) are not essential 
as long as the imaginary part is finite. 

The counter example above is not only of academic interest. Schweitzer and CzychoU observed as well that the 
inclusion of a nearest-neighbor self-energy leads to wrong signs of the imaginary parts. They considered the 1/d 
expanajpD of a C/^ perturbation theory around Hartree-Fock for the Hubbard model and the p£iriodic Anderson 
modeB'tZl. They reached consistency by including higher l/d corrections (for d = 1 up to 50 terms )la. Problems with 
the aiialyticity (uniqueness) of the solution occurred also in the first investigations of 1/d corrections in the Hubbard 
modeH (Falicov-Kimball modelEfl). 

To the author's knowledge there is no necessary or sufficient theory so far, which predicts under which circumstances 
such problems have to be expected or can be excluded. A sufficient argument excluding wrong signs of the imaginary 
part of the approximate self-energy is given by the theorem: 

If the approximation considered can be interpreted as an expansion of the self-energy in a parameter A 
and if m is the leading order, in which the imaginary part of the self-energy does not vanish, then the self-energy 
approximated in the m-th order has the right sign. 

The proof relies on the continuity of limits if the expansion exists. According to the precondition holds 

> ImEA(w,k) = A'"ImS(™)(w,k) + 0(A("+i)) , (3) 

which is equivalent to 

0> lim A-"TmI]A(t^,k) = ImS(™)(tj,k) . (4) 

The index k is the wave vector in a homogeneous, translationally invariant phase. The derivation for general phases, 
for instance the AB-CDW, is given in appendix A. 

The derivation in (^) and in (^) holds strictly only for the non self-consistent treatment. In the generic situation, 
however, the leading order of the self-energy with non-vanishing imaginary part results from a certain diagram class 
and the analytic properties do not depend on the specific form of the Green function entering. If this is the case, 
the statement of the theorem extends also to the self-consistent treatment where the quantitative form of the Green 
functions are not known a priori. 

The theorem helps one to understand the observations made by Schweitzer and CzychoU. In the expansion 
of the d-dimensional Hubbard model and of the periodic Anderson model one has A = 1/rf and to = since the 
self-energy is imaginary already in the first order. For the perturbation theory in V one has \ — U and m — 1 
since the self-energy stays real in Hartree-Fock. Applying the rationale of the theorem twice one understands that 
the self-energy in U"^ of the infinite dimensional model has the right analytic behavior. If further 1/d corrections are 
included this does not need to be true. The result of Schweitzer and CzychoU, that the linear 1/d correction leads 
to wrong signs, proves that the theorem is sharp: If the precondition fails, the implication fails, too. The second 
obervation, that the inclusion of very many 1/d correction terms remedies the failure, can also be understood easily. 
In this case the calculations approximates the U'^ perturbation theory of the finite dimensional models very well. 
According to the theorem, this perturbation theory displays the right sign, too. 

The_above observations indicate that also the analyticity problems encountered for 1/d corrections in the Hubbard 
modeld are not due to the approximations used to solve the effective impurity problems. Rather each time that the 
theorem does not apply one has to expect that analyticity problems arise for certain parameters. Considering eq. (6a) 
in ref. 18 or equivalently eq. (370) in ref. 3 one realizes that the spectral density of the local self energy might change 
sign. This cannot obviously be excluded from the way how the impurity self energies are computed. 

Turning to the 1/d expansion of the present model of spinless fermions (A = ^/d), one notes that the theorem 
applies with to = 1. Therefore, the equations including linear 1/d corrections display the right analyticity. These 
equation will be set up in the following. 



A. Resulting equations and one-particle results 

This section is kept very concise since it contains material which is partly published elsewherei. For two reasons, 
however, it cannot be omitted. Firstly, a different notation using different intermediate quantities shall be introduced. 
Secondly, the one-particle results are necessary requisites to understand the conductivity results in the subsequent 
section. 

The treatment of a self-energy of the type depicted in fig. || is commonly known (see e. g. refs. 11, 19, 20). Dealing 
with the symmetry broken phase, however, requires some extension. In a previous worm local Green function and 
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the self-energy are distinguished according to the sublattice to which they belong. In the present work, sum and 
difference of the quantities on the two sublattices will be used. The local quantities on site i belonging to sublattice 
T e {A, B} are 

gr := Gi^i{uj) (5a) 
S.H:=SH(^) + SC(,;) , (5b) 

where Gi^i{uj) is the full local Green function and S is the local self-energy. The Fock part will be treated subsequently. 
The index ^ stands for the Hartree term (first diagram in fig. ||); the index ^ stands for the local correlation (third 
diagram in fig. ^). Let us define 



5s('^) 
S(^) 



= igAito)+gBiLu))/2 (6a) 

= {gA{u;)-gB{io))/2 (6b) 

= (SA(w)-f Sb(c^))/2 (6c) 

= (Ea(co) - I]b(c^))/2 . (6d) 



The spectral functions of the Green function are called A^s and N^y, respectively; the spectral functions of the self- 
energy S and A are called iVs and TVa, respectively. The non local Fock term is := j-, where i and j are adjacent 
sites. It turns out, that is negative (for repulsive interaction), real, and that it does not depend on whether the 
fermion hops from A to _B or vice versa. Hence, it renormalizes the hopping 

t ^ jt with 7 := 1 - \/ZJ:^/t . (7) 

Note that for attractive interaction 7 could become which would lead to a breakdown of the theory. Such a singularity 
is absent in the repulsive case. 

In the AB-CDW, the modes at k couple to those at k + Q. Hence one has 



Gk+Q,k Gk+Q,k+Q / I -A(tj) CJ-I](cj)+7£ 



(8) 



From this equation one obtains 



5s (^) = — / 9^^,/ ; go(\/w^ - A^(a;)/7) 
^\Jw^ — A^(a)) 



(7^)2 - A2 



No{e)d£ (9a) 



A(^) 



^D^'^) = / / A 9/ ,- go(V^''- A2(a;)/7) 
^\Jw^ — A-^(Ci;) 



,2 .--^^2 ^^o(e)d£, (9b) 



lip- — (72)2 — A2 



-C30 



where w is short hand for — E(aj). l— , 

The averaged Hartree term U{nA + 't-b)/2 renormalizes the chemical potential^. The Hartree contribution to 

A is Ub where b := (riB — np^jl is the order parameter, i.e. the particle density difference. It is given by & = 

00 

— J NY)(uj)f-p{uj)dLu, where /f(w) is the Fermi function. 

—00 

The Fock term can be calculated from the nearest-neighbor Green function Gj+oj 

oc 

= / lm{G,+a,j{uJ + 0i))fF{io)dLO, (10) 



which is given by 



G,+a,(-) = --^|£(k)Gk,k(gp 

z 

[{Lu~i:)gs{Lo)-Agu{Lo)]ds. (11) 
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The Fock term is related to the kinetic energy = {U /Z^/'^){T). Thus, (10) can be evaluated using ( pi] ) and (||). 
The local correlation term is given in terms of the Matsubara frequencies ljJ\ (fermionic) and lji (bosonic) by 



(12) 



Here, the index r stands for the other sublattice, i.e. for ^ if r = i? and vice versa. By performing the Matsubara 
sum one obtains the convolution 



oo oo 



N-{uj")N-{uj" - uj')Nr{u: - uj') 



— OO —OO 



- u)h{-uj")h{u:" - u') + /f(co - uo')h{u:")Mu:' - uj")] duo' dw" 



(13) 



for the spectral function N^^{llj) belonging to E^(aj). The convolution can be expressed most conveniently in the 
Fourier transforms 



N'^it) := / exp(-iwt)iV(±w)/F(-w) 



N{t) := / exp {-iut)N{uj) . 



(14a) 
(14b) 



Eq. ( p^ ) becomes as simple as N^^ (t) — ^ 



N+N+NZ 



+ N-N-N+ 

T T r 



. In sums and differences one obtains 



NAit) 



IP r 

Z 



C/2 r 



(15a) 
(15b) 



The complete self-energy E and A are given by the following inverse Fourier transforms 



T,{lu + Oi) ^ exp {iujt - Ot)N^{t)dt 



A{lu + 0i) = Ub- i J exp {iujt ~ Ot)NA{t)dt . 





(16a) 
(16b) 



In (16b) the Hartree part has been added. 

So far, no assumptions concerning the DOS entered. The formulae hold for all fillings. At the particular value of half- 
filling the additional symmetries Ns{uj) = Ns{—uj), A^d(^) — —Nui—uj), N^{lu) = N^{—u!) and Na{uj) = —Na{—co) 
can be exploited. The fact that the spectral densities are real tells us that N{—t) is the complex conjugate (c.c.) of 
N{t). Thus (ph simphfies at half-filhng to 



iVs(i) = 

NA{t) = 



El 
Z 

El 
z 



{iKr-{Nir}N+ 



c.c. 



(17a) 
(17b) 



This terminates the set up of the equations which have to be solved self-consistently on the one-particle level. 

For those who intend to implement these equations or similar ones some remarks on the numerical realization are 
in order. As usual, the self-consistent set of equations is solved by iteration. At T = it is favorable to use a relaxed 
iteration. This means that the self-energy E and A from the n-th and from the n + 1 iteration are averaged and used 
for the subsequent calculation instead of using only the n -I- 1 iteration. This procedure damps oscillatory deviations 
from the fixed point more rapidly. It is even more advantageous to let the programme decide whether relaxed or non 
relaxed iteration converges faster. 
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The Fourier transformation is the most time consuming step. The best algorithm for this task is the so called 
Fast Fourier Transformation (FFT). The extremely large number of points, which can be used with the FFT, over- 
compensates the disadvantage of an equidistant mesh which cannot be adapted to regions where the DOS changes 
rapidlju. In the AB-CDW 2^^ points were used. The vectorization on a IBM3090 still permitted to do one iteration 
step comprising four FFT in 19 seconds. A very good precision could be achieved. The sum rules 



r/^ 1 /I 



oo 

r rP /I 

J NA{L0)dL0^—bi^--b' 



(18a) 
(18b) 



are preserved up to 10" . Note that (18b) holds only at T = whereas (18a) holds for all temperatures. 

In order to achieve the high precision also at T = 0, it is necessary to discretize the DOS carefully. At the gap edges 
the DOS displays inverse square root divergences a/y/to — wa- The parameters a and uja are determined directly 
from the self-energy using (^. The diverging part of the DOS is discretized by using the average value in the interval 
[uji — Su!/2,uji + S(jj/2] instead of the DOS value at tOi. 

Once the Fourier transforms are essentially linear one as to avoid a non-linear time loss in the calculation of the 
complex free Green function go{z). Therefore, the integration from the Hilbert representation must be avoided. This 
is done by using the approximate expression 



r 13033 
\ 29088 



8675 



174528 

r4167 459 729 ^ 

\6464 6464^6^ ^ 12928^ 
4167 459 729 ^ 

6464 ~ 646476^ 12928^ 



v/6 - £2 - 



2/3- (e- 2^/2/3)2 - 
|\/2/3-(e + 2v/273)2 



(19) 



for the three dimensional DOS N3{e). The identities h{z;a) := (I/tt) Va — £^/(^ ~ £)d£ = z ± \fz^^~~a and 

(I/tt) \'f^^E^f^i~^l{z - e)de = -(a/2) + zh{ z;a) permit to compute go{z) for any z quickly. The r.h.s. of ( p^ ) is 
chosen such that the van-Hove-singularities are at the right places and such that the first moments (including the 
8th) are reproduced exactly. The relative accuracy achieved is 4 • 10~^ for A'3(0) and 10~^ for the 10th and the 12th 
moment. 

The calculation of the Hartree and of the Fock parts are linear in the number of discretization points. Concluding 
the remarks on the numerical realization we state that all parts of an iteration step are essentially linear in the number 
of points used. This allows a reliable and efhcient computation. 
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FIG. 6. Density of states and spectral function of the self-energy in the homogeneous phase at f7 = 2 and T = and T = 2 
in d = 3. For definitions see eqs. (n). 
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In fig. |6|, results for the DOS and the spectral density of the self-energy in the homogeneous phase are shown. The 
spontaneous symmetry breaking is deliberately suppressed. Only positive frequencies are displayed since the functions 
are even. At T = 0, one notes that the imaginary part of the self-energy tends quadratically to zero for w — > 0. From 
( |l3| ) this follows for all free DOSes with finite non-singular value at the Fermi edge. Tlaus the homogeneous low 
temperature phase of interacting spinless fcrmions is a Fermi liquid. But this phase is thcrmodynamically unstable 
(see below). The DOS still bears signs of the van-Hove-singularities which are smeared out only a little due to the 
interaction. Note that the width is increased by the Fock term. In the free case the half- width is V6 « 2.45. High 
temperatures smear out the minimum of iVs at a; = completely. The solution depicted is stable since at T = 2 no 
AB-CDW is possible. 
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FIG. 7. Density of states and spectral function of the self-energy in the AB-charge density wave at C/ = 2 and T = 
(6 = 0.311005), T = 0.225658 (6 = 0.250000) in d = 3. For definitions see eqs. (|). The sum quantities in (a) and (b) are even 
functions of frequency; the difference quantities in (c) and (d) are odd functions. 



In fig. 0, stable solutions with b > are shown. Note the square root divergence in the DOSes (left column) in 
the vicinity of the gap. At T = the gap is at 2u!a ~ 0.6 whereas the spectral density of the self-energy becomes 
finite at about 1.8 ~ 6uja- This results from the two convolutions involvedB. They make the gap in the density of the 
self-energy to be exactly three times the gap in the DOS. Put differently, the finite spectral density of the self-energy 
corresponds to the inelastic scattering of a particle or a hole involving an additional particle-hole pair. Thus, the 
necessary minimum energy is three times the elementary gap. The physically important implication is the existence 
of quasi-particles with energies between and Swa with infinite life-time. Following the arguments of LuttingertiJ 
by which he shows that the density of the self-energy generically goes like cj^ at the Fermi edge one comes to the 
conclusion that this factor 3 is not an artifact of the approximation but valid to all orders. Therefore, if the conditions 
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are such that the the homogeneous phase is a Fermi Hquid, i.e. Luttinger's argument holds, a gapped, spontaneously 
symmetry broken phase has a factor 3 between the gap in the DOS and the gap in the self-energy. This implies 
also the existence of undamped quasi-particles which have interesting consequences on the transport properties (see 
below). The exponent of the power law with which the imaginary parts of the self-energy rises at w = Swa is 3/2. 

At finite temperatures the energy gap is smaller since the order parameter has decreased. This effect is visible 
already in the Hartree treatment. In addition, the energy gap is smeared out: thermal fluctuations represented by the 
local correlation term S'-' induce a certain spectral weight within the "gap" which docs no longer exist in the rigorous 
sense. The occurrence of two maxima in Ny: and in should be noted. 
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FIG. 8. Density of states and spectral function of the self-energy in the AB-charge density wave at f/ = 8 and T = 
(6 = 0.479312), T = 1.509384 (6 = 0.260004) in d = 3. For definitions see eqs. (||). At tj « 12.0 hardly visible satellite bands 
are present in A^'s and Ad for T — 0. They result from the imaginary parts of the self-energy around this frequency. 

In fig. H the generic results for large values of the interaction are shown. At T = the factor 3 between the gap 
in the DOSes and the gap of the spectral densities of the self-energies is even more easily discernible. At the finite 
temperature (T « 1.5), all the structures are smeared out; the order parameter is considerably smaller than at T = 0: 
b = 0.260 at finite T to 6 = 0.479 at T = 0. The comparison of the spectral weights of the self-energy at zero and 
at finite temperature illustrates an important effect. The correlation term is suppressed by the symmetry breaking. 
The larger b the smaller is the area under the curves in fig. |^(b) and (d). The effect can be understood quantitatively 
with the help of the equations ( [l8| ) which imply that the area under the curves vanishes for 6 ^ 1/2. This leads to 
the counter-intuitive effect that the significance of the correlation term decreases on increasing interaction at T = 
albeit it is quadratic in the interaction 

In fig. H, hardly discernible satellite bands exist at u; « 12. They are engendered by the finite imaginary part of 
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the self-energy at these energies (see fig. ||(b) and (d)). To demonstrate that there are in fact infinitely many satelhte 
bands with exponentially decreasing weights, the densities and N^: are plotted logarithmically in fig. ^. 
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FIG. 9. Density of states A^s (short dashed curve) and spectral function A^e (long dashed curve) in the AB-CDW at f7 = 8 in 
logarithmic scale. The difference quantities are not shown since their values lie only slightly under those of the sum quantities. 

The principal band of the DOS consists of quasi-particles with infinite life-time at loa ~ 4. The satellite bands 
correspond to peaks in the spectral density of the self-energy. The satellite bands are located at (2m -I- l)a;A where 
m is an integer. The peaks in the spectral density of the self-energy are located at (2m + 1)u!a where m is an integer 
but not or —1. This phenomenon is generic for the self-consistent solution of a system of equation comprising 
convolutions of strongly peaked functions. It appears only at large values of U because it is necessary that wa ~ U/2 
is larger than the band width in order to resolve the peaks. Note that according to (^), a large value of A induces 
band narrowing. Whereas the principal band is y/Q wide at J7 = 0, its width is shrunk to about unity in fig. ^(a). 

For detailed numerical results on the order parameter as function of interaction and of temperature as well as on 
the critical teirmepature the reader is referred to ref. 6. The asymptotic behavior at small U is discussed analytically 
by van DongenE3^c3. In a nutshell, the correlation term renormalizes the Hartree results for b and Tc by a constant 
factor of order unity which tends to unity for d oo. 



III. CONDUCTIVITY: FOUNDATIONS 



Due to the point symmetry group of the hypercubic lattices the conductivity cr(w) can be treated as a scalar. 
Previoius one-particle results showed that the treatment on the level of linear 1/d corrections should yield reasonable 
resultgj in d = 3. 

The conductivity is calculated from a two-particle correlation function. This will be done here from the current- 
current correlation function x"^^ ■ The conductivity comprises two contributions (t{lu) — (Ji{u!) -f (T2('-^)- The ficst term 
depends on the occupation of the momentum states (nt) whereas the second term is proportional to x''''(cli)E3 



-iM--y (20a) 

BZ 

^2{to) - -x"{^) . (20b) 



The current-currrent corijelation function will be computed including 1/d corrections with the help of the 
Baym/Kadanoff formalisnii3'Ej. Specific correlation functions are determined from the general two-particle corre- 
lation function L(12, 1'2') via 



= J A{l,l')L{l2,l'2')B{2,2')dn'22' . (21) 



The numbers stand for composite space and time coordinates (or momentum and frequency coordinates). The 
measure dW 22' tells which coordinates are integrated. The quantities A and B represent the operators for which the 
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FIG. 12. Current-currrent correlation with 1/d corrections. 
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caused by the non local Fock term is left. For comparison: in the infinite dimensional Hubbard model the simplifi- 
cations are even more drastic. All vertex corrections drop out and the current-current correlation function is just the 
convolution of two Green functionsE^. 

Let us call the value of the first diagram in fig. [ij the "free" current-current correlation function and let us use the 
symbol Xo^ for it- In the homogeneous phase one obtains with the help of (^l|), (24) and of the propagator in k-space 
(c. ~ EH - 7£(k))~i 

jj _ ^ _ 4T ^ f sin2(fci) 



LUi, —UJx—UJ„ 



^ ,„ J^-, , J (w^- 7e(k)) [wx - 7e(k)) {2Trf 



BZ 



where w^/x := iw^/A - 

We focus now on the segments between two wavy lines in fig. The conservation of energy and of momentum 
makes it possible to carry out the sum over all momentums and energies by considering independent momentums 
and energies circulating in each segment. Then the momentum in a wavy line is the difference of two adjacent wave 
vectors k and k'. A second time, the evenness and the oddness in the components of the wave vector is used to write 
for the factor of an interaction line 

U ^ 2U v-^ 

- — ^ cos(fcj - kl) = — [sin(fcj) sin(A;-) + cos(fci) cos(fc-)] 

i=l i=l 

2U 

sin(A:i)sin(fci) . (26) 

The argument is obvious for one of the border segments and follows for those in the middle by induction. 

At the end one realizes that each segment corresponds to a factor of —{U/2)xo'^ which justifies to call the right side 
of fig. O a geometric series which takes the value 



X (w + 0^ = jj, — , ,^ (27) 
1 + Uxo + 0i)/2 



after analytic continuation. The derivation of a similar formulae in the AB-CDW is given in appendix B. The results 
are cited below. 

The momentum integration in ( p5| ) requires a modified DOS, to be called the conductivity DOS henceforth 

sin'{k,)S{u;-s{k))^ , (28) 

BZ 

oo 

from which we define also the conductivity Green function gc,oiz) J Nc,q{uj) I {z — Lu)duj. The conductivity DOS 

— oo 

can be simply derived once the DOS is known. These two functions are related via 

Noicj) = --^(c.) . (29) 
uj ouj 



This relation stems from the fact that one has to replace one of the d factors (l/7r)l/\/t^ — uj'^ in the convolution 
for the DOS by (l/7r)vi^ — in order to calculate the conductivity DOS. The derivation uses the representation of 
convolutions as products in Fourier space. 

Using the definition of the conductivity Green function and partial fraction expansion it is straightforward to rewrite 

Xo («^m) = ^ > - (30) 

7Z Wi, ~ Wx 



Analytic continuation of the latter gives the general formula (eq. (14) in ref. 26) for the current-current correlation 
function in the homogeneous phase. 

In the AB-CDW, it is also possible to sum the series in fig. |l^ as geometric series. The main difcrence is the fact 
that 2x2 matrices instead of scalars are involved. The details are given in appendix B; the resultsEZl are 
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where the quantities Ai , A2 and A3 are defined by 
2UT 



Ai{iu}„ 



E 



A2{iuJm) 



2UT 



E 



(WA - W^){gc,s{'i^iy) + .gc,s(«WA)) - (A(zWa) - A(iWj.))(gc,D(«Wi.) + 3c,D(iWA)) 



A3 {iuJn 



2UT 



E 



A(wA)gc,s(^^'>) ^ ^xgcpjii^i.) + A(^u;^,)ffe^s(^tx^A) - Wivgc,D(»^A) 

wl - wl - iA^iu,) - A^iux)) 

In complete analogy to the usual Green functions, the conductivity Green functions are gc,s ■= (<?c,A 
ffc,D := (ffc.A - 5c,b)/2, hence 



5c,s(t^) 



^\Jw'^ — A^[ijj) 

which compares to (^) (w is short-hand for lo — E(cl;)). 

Now a relation for the dc-conductivity shaU be derived. In order that the limit hmi^^o exists 



(32a) 

(32b) 

(32c) 
gc,B)/2 and 

(33a) 
(33b) 



^ (O^-iaAf (2^ - V 

BZ 



(34) 



must hold according to (|2y). The operator T stands for the kinetic energy. Eq. (^4|) implies also the /-sum rule 

00 

/ (ix"'"'/^)'^'^ = —■n{T)/d. At the end of appendix B, it is shown explicitly that ( p4| ) is valid since A3 vanishes at 

—00 

LO = Q and Ai = — [/ (T) / (270?) = 1 — I/7. For the dc-conductivity one obtains 



a(0) 



dx 



jj 



2^7^ aAi 



cJ = 



(35) 



j=0 



For explicit evaluation it is useful to split tTdc(O) into a term including retarded and advanced Green functions (7dci 
and a term including only retarded or advanced Green functions (Tdc2 after analytic continuation. This yields 



Cdcl 



ttZ 



(1 - Rc^)N,^s - (ReA)jV,,D 
(1 - ReE)Afs + (ReA)AfA 



{-Uuj))duo 



(36) 



where f^i^) is the derivative of the Fermi distribution, and 



0"dc2 



7 



(1 - T,)d^gcs - Adujgc.B 



ui+Oi 



1 



(1-S)(9„E-1) + A9<,A 
1-Re / ((c. - I])5s - A5d)<,+o, (-/fM)^;'^ 



(37) 



In the last expressions, all the Green functions are retarded. In the homogeneous phase, the contribution ( |3^ ) is more 
important than the one in (^7|). The former diverges for T and uj —> 0, the latter does not. In the symmetry 
broken AB-CDW, however, both terms turn out to be essential. 

Eqs. (|3l|), (^), ( ^ ) and ( |3^ ) are the foundation for the calculation of the conductivity for zero and for non-zero order 
parameter. The focus of the present work is on the AB-CDW. The properties of the conductivity in the homogeneous 
phase (e.g. Fermi liquid behavior) are presented in detail in ref. 26 where also the influence of the truncation of the 
1/d expansion is discussed. 
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IV. CONDUCTIVITY: RESULTS 



In this section we present and discuss results which follow from the general equations derived in the previous section. 
All results are calculated at half-filling and for d = 3. 




FIG. 13. Scaled real part of the dynamic conductivity Re (7{u)) in the non symmetry broken phase ai U — 4.243 for T — 0.393 
(solid lines), T — 0.196 (dashed lines), T — 0.049 (dashed-dotted lines), and T = 0.025 (dotted lines). Main figure: spinless 
fermions in d = 3; inset: Hubbard model in the non-crossing approximation (data from Th. Pruschke). 

In fig. |l^, the real part of the dynamic conductivity is depicted in the non symmetry broken phase for different 
temperatures, i.e. the occurrence of a symmetry broken phase at low tertmerature is discarded deliberately for the 
moment. They are compared with results of Pruschke, Cox, and JarrellESO for the half-filled Hubbard model in 
d = oo, obtained in the non-crossing approximation. In both cases the interaction value is ^-t= 4.243 (in our units) 
which is just below the value where the Mott-Hubbard transition occurs in the Hubbard modelEj. For spinless fermions 
the Drude peak is absolutely dominant. Its weight is very large. Its width is given by the imaginary part of the self- 
energy at the Fermi level N^{0) (see ( [3^ ) with A = or eq. (14) in ref. 26) , i.e. the width is proportional to T^. The 
shape of the Drude peak corresponds very well to a lorenzian. 

Only at low temperatures a shoulder emerges. This shoulder is the effect of interaction induced scattering. The 
fluctuations are not particularly strong. It was already shown previouslycJ that the average over the Z interaction 
partners reduces the relative fluctuations. There is no Mott-Hubbard transition without symmetry breaking in the 
spinless fermion model because an increasing interaction enhances not only the fluctuations but also the Fock term 
(absent in the Hubbard model) which stabilizes the Fermi liquid phase. These features are particularly obvious in the 
comparison with the Hubbard model data. In this model, the Drude peak is very reduced at all displayed temperatures 
since much of the weight is shifted to the peaks induced by the strong local particle density fluctuations. 

Besides the difference shoulder vs. peak it is interesting to note the difference in energy scales. In the Hubbard 
model, it is more or less U which sets the energy at which the peak occurs. This can be understood as the energetic 
effect of whether or not an electron with a different spin is present. The typical energy for the shoulder is obviously 
much smaller. This in turn can be understood in the same way as before but it has to be taken into account that the 
number of possible interaction partners Z leads to a reduction of the relative fluctuations of the order of 1/ \/Z. This 
yields an energy of roughly 1.7 in the particular example which is in good agreement with the numerical result. 

Due to the nesting at half-filling, the system of spinless fermions undergoes a transition to a spontaneously broken 
translation symmetry for all (positive) values of the interaction on lowering the temperature. This spontaneously 
broken discrete symmetry implies the occurrence of a gap which grows exponentially wa c>c exp(— c/Ii.) for low values 
of the interaction at T = (see ref. 6 and refs. therein). It is visible in the dynamic conductivitjitJ. In fig. |lj, its 
growth on decreasing temperature is shown in four snap-shots. 



15 





FIG. 14. Real part of the dynamic conductivity Re o-{lj) in ci = 3 at 17 = 2.0 in logarithmic scale. Fig. (a) T — 0.300000 and 
6 = 0; fig. (b) T = 0.225658 and b = 0.250000; fig. (c) T = 0.155286 and b = 0.299801; fig. (d) T = and 6 = 0.311005. In fig. 
(d) the (5-peak at uj = 0.13312 is not shown, its weight is 0.062336. 

In fig. |lj(a), T is still above its critical value. No structure is visible except for the dominant Drude peak already 
discussed in fig. |l3|. In figs, ^(b)-(d) the gap is present and discernible. Its value is approximately 2wa if i^A is the 
value of the energy gap in the DOS, see figs. and ||. But there is also some weight within the gap for T > since 
the correlation contribution blurred already the gap in the DOS. Note in passing that the f-sum rule can be verified 



numerically on the results shown in fig. 14 very accurately (to the fraction of a percent at T = 0; to the fraction of a 
permille in the homogeneous phase). 

The Drude peak does not vanish immediately in the AB-CDW. It becomes smaller and narrower on decreasing 
temperature. Its maximum value does not vanish for T — » (see below) but its weight does. In fig. nS, 
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FIG. 15. Enlargements of two frequency intervals for T — 0.0833833 and 6 = 0.310773. Fig. (a) shows details of the Drude 
peak; fig. (b) the excitonic resonance. 



two frequency intervals are shown in detail for a fairly low temperature. Fig. p^a) displays the Drude peak again. 
The interesting feature is its small width (compared with the width of the Drude peaks in figs. |lj(b) and (c)). It 
cannot be explained by a factor of but corresponds to an exponential shrinking exp{—uj/\/T). As already observed 
in the one-particle properties, an increasing gap reduces the influence of the fluctuations. 

Fig. |l5|(b) shows a very interesting feature below the proper band edge at w « 2uja- This resonance is also visible in 
fig. |lj(c) whereas the resonance and the band edge are not resolved at a higher temperature, fig. |l^(b). The resonance 
can very well be approximated by a lorenzian. At T = 0, it is also present as a (5-peak (not shown in fig. ^(d)). It 
originates from a zero of the denominator in (^). At T > 0, only the real part of the denominator vanishes and its 
imaginary part leads to the observed broadening which depends strongly, namely exponentially, on the temperature. 

Physically the resonance can be interpreted as a bound state, an exciton, between a particle in the upper band and 
a hole in the lower band in the reduced Brillouin zone of the AB-CDW. The energy difference between the position 
of the exciton and the band edge is its binding energy. The type of diagrams which yield the denominator in ( |3l| ) 
corroborates the interpretation as an exciton. The vertical interaction lines stand for the repeated interaction between 
particle and hole in the two propagators involved in the calculation of x'^'^ ■ It should be noted that, for instance, for 
the parameters of fig. ^(d) about 70% of the weight of the conductivity are found in the excitonic resonance (one 
may not be misled by the logarithmic scale). This means that the excitonic effect is not at all a small side effect. 

Concluding the part on the dynamic conductivity, we discuss fig. |l^ 
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FIG. 16. Real part of the dynamic conductivity Re i7{uj) in d = 3 at (7 — 8.0 for T = in logarithmic scale. The (5-distribution 
is not displayed. 

which shows results for a large interaction value U. Due to the induced large gap and due to the narrow effective 
band width several frequency intervals of absorption are well separated. The peaks are caused by the convolution of 
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the satellite band presented for the one-particle properties. Note, however, that the weight of these satellites decreases 
rapidly by a factor of 100 from peak to peak. These small amplitudes render an experimental verification certainly 
extremely difficult if not impossible. Nevertheless, it would be interesting to know whether such satellites exist. Their 
existence would support the application of a self-consistent approximation since the non self-consistent calculation 
yields only two peaks besides the (5-peak which is not shown. 

Since the dc-conductivity in absence of symmetry breaking has been extensively discussed in ref. 26 we will treat 
here exclusively the case with symmetry breaking. The result of (^) and (|3^) is depicted in fig. ^ 
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FIG. 17. Temperature dependence of the dc-conductivity at U — 1.0 (fig. (a)) and at ?7 = 8.0 (fig. (b)). Below T = 0.026 in 
fig. (a) and below T = 0.6 in fig. (b) a fit was used (see main text). 

for weak and strong interactiorH. Tp_the right of the cusp the system is in the non symmetry broken phase. The 
conductivity is essentially proportionalEj to T^. On entering the symmetry broken phase with gap, the conductivity 
falls drastically since the energy gap reduces the DOS at the Fermi level. Surprisingly, however, the conductivity 
does not vanish for T ^ although the DOS vanishes in this limit. There is even a very slight uprise of adc close to 
T = 0. This phenomenon is again a manifestation of the suppression of correlation effects by the energy gap. The 
DOS is reduced by a factor of exp(— ojA/r) but so is the imaginary part of the self-energy in (32) which is responsible 
for the quasi-particle life-time. These two effects cancel exactly. Put differently, an exponentially small number of 
quasi-particles of exponentially large life-time carries a constant current (but see discussion below). It remains an 
algebraic dependence on T of the dc-conductivity The constant term and the linear one can be computed analytically 
and where used to complete the curves in fig. O for small values of T where the numerical calculation is no longer 
precise enough due to extinction. 

The limit value limT^o <^{^ — 0) is given in fig. |l^ 
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FIG. 18. dc-conductivity a{u! = 0) in the limit T ^ in logarithmic scale. 
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as function of U. As expected it decreases rapidly for U ^ oo. Note the logarithmic scale. 

What do the above findings for a^c imply for the existence of a metal-insulator transition? Seemingly, even 
spontaneous symmetry breaking does not suffice to render the system insulating. But it must be noted that the 
"residual" conductivity limr^o "'(^ — 0) is infinitely fragile: any other arbitrarily weak scattering mechanism which 
does not die out on T — > e.g. disorder or scattering at the borders of the sample, will take over. The exponentially 
vanishing DOS will yield an exponentially vanishing dc-conductivity. This is reflected in the exponentially decreasing 
width of the Drude peak which, at constant height, implies an exponentially decreasing weight. 

Experimentally, very pure samples might allow to see the beginning of the plateaus in fig. |l^ before the above 
cited other scattering mechanism reduce the conductivibi^ This behavior is in complete analogy to the one observed 
for the shear viscosity ri{T) of Helium 3 in the B phasetj. In this system like in the system of spinless fermions in 
the AB-CDW one observes an exponentially diverging mean free path since the collision between (quasi-)particles is 
suppressed by a gap. In the so-called "Knudsen regime" collisions of quasi-particles with the wall of the container 
dominate the collisions between the quasi-particles. In Helium 3, one observes a sharp drop below Tc and then the 
beginning of a plateau before finally r]{T) vanishes rapidly. The theoretical result for the infinite system predicts a 
gentle uprise just like the one we prec^t in fig. In both cases, a factor exp(— wa/^) in the DOS cancels with 
the same factor in the scattering ratec2l. This interesting analogy underlines the validity of the results of our 1/d 
approach. 

V. DISCUSSION 

Two main questions arc addressed in the present paper: (i) How one can an infinite dimensional result be improved 
by including 1 /d corrections in a systematic way? (ii) Which influence does spontaneous symmetry breaking have on 
the conductivity? 

It turned out that it is highly non trivial to construct systematic and reasonable approximations to arbitrary 
order. This is true already on the conceptual level. It was argued in detail that the self-consistent calculation has 
certain advantages since it yields thermodynamically consistent and conserving approximations. The Baym/Kadanoff 
formalism, however, is not sufficient to guarantee an approximation which is free from obvious contradictions. It 
was shown that an inappropriate approximation may lead to the wrong analytic behavior of Green functions and 
self-energies even though the approximation was derived from a generating functional. 

A general theorem was presented which allows to judge whether wrong analyticity may occur. If the conditions 
of the theorem are fulfilled the appearance of the wrong analyticity is excluded. This theorem explains a couple of 
observations which were made in the last years on the application of perturbation expansions and/or \/d expansions. 
It is used to show that the self-consistent treatment of 1 /d corrections for spinless fermions is a good approximation: 
it possesses the necessary analytic behavior. 

For 1/d corrections in the Hubbard model the presented theorem makes no statement since the self-energy has 
already an imaginary part for d — oo. This does not imply that the systematic inclusion of 1/d corrections for the 
Hubbard model is impossible, but one may expect further difficulties. As a matter_pf fact, analyticity problems have 
been encountered in the first calculations of 1/d corrections in the Hubbard modelB. 

It should be stated that the self-consistent treatment of l/cLjcorrections to any finite order in 1/d remains a mean- 
field theory. As in the d = oo treatment of the Hubbard modelEI the mean field is dynamic, i.e. it retains a dependence 
on frequency. But in the skeleton diagrams, which are considered in any finite order in l/d, only lattice sites of finite 
distance occur. This means that critical fluctuation&-|are always cut off. In d = 1, for instance, the inclusion of 1/d 
corrections reduces the order parameter considerablyB but does not destroy the order completely. 

In the self-consistent 1/d treatment of spinless fermions, two-particle properties can be reached, too. In this work, 
the Bethe-Salpeter equation was set up in general and solved in the particular case of the conductivity cr(w). This was 
possible for the non symmetry broken phase as well as for the charge density wave. The equations were evaluated in 
d = 3 since the approximation should yield the best results for this value of all experimentally accessible dimensionsO. 

A number of phenomena were described in the 1/d expansion which can be compared with other theoretical 
predictions or experiments: 

- the dynamic conductivity cr(w) in the homogeneous phasc-bas a Drude peak. Its width decreases quadratically in 
T for small values of T. The dc-conductivity is always finiteEj. 

- The Drude peak persists in the CDW but its weight vanishes exponentially cx exp{—Lu^/T), where loa is the gap in 
the one-particle spectra. The height of the Drude peak, however, does not vanish since the diverging quasi-particle 
life-time cancels the vanishing density of states. 

- The real part of cr{uj) displays a band edge at w 2uja- The singularity at the edge is a square root. Just below the 
edge an excitonic resonance is situated which is the bound state between a particle and a hole in the empty and in 
the full band, respectively. These bands are created by the spontaneous symmetry breaking. 

- For strong interactions the real part of (T(a;) shows exponentially decreasing peaks at cj 2mwA; w G {1, 2,3,.. .}, 
which reflect the peaks in the one-particle DOS at « (2m — 1)u!a- 
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- Strictly speaking, there is no metal- insulator transition. But the Drude weight decays rapidly on T ^ 0. Finally, 
other scattering mechanisms will dominate over quasi-particle-quasi-particle collisions. 

In summary, we conclude that the self-consistent treatment oil/d corrections describes successfully a large variety 
of phenomena since it includes the leading frequency dependence of the self-energy. It is a generalized and improved 
mean-field theory. 
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APPENDIX A: GENERAL DERIVATION OF THE THEOREM 



Let us consider a general correlated fermion problem without magnetic field. A one-particle basis {h^} {bf fermionic 
creation operator) can be chosen in which the free (one-particle) Hamiltonian can be represented as a real matrix 
h^^'' . The retarded self-energy and the retarded full (interacting) Green function are matrices as well. They might be 
complex. The argument runs at an arbitrary but fixed value of uj. According to the Dyson equation one has 

.(0) 



Lo + Oi- h^"> - Sj . (Al) 

The full Green function can be written as G = Re G + i Im G, where the real and the imaginary parts are real, 
symmetric matrices. The same is true for the self-energy. Let us define in particular B .= ImS. 

The aim is to show (a) that the leading order B*-™-* in A of the matrix B is positive semi-definite and (b) that this 
implies that all spectral densities are non negative. The main difference to the argument in the main part is that the 
matrices do not commute in general. I— . 

All one-particle Green functions — < T {a{t)a^ {^)} > have positive spectral densitiesEj. Here T is the time ordering 
operator and a"*" (a) is an arbitrary fermionic creation (annihilation) operator. Hence, one has for the corresponding 
vector V defined by a"*" =: Vibf 

< v+ (Im G) V . (A2) 
Since the above equation holds for any v, it implies that the imaginary part of G is negative semi-definite. Expanding 



( Al ) in powers of B and resumming the imaginary part yields 

ImG = -A-^B A-^l + B A-^B A-^) , (A3) 

where A := ui — h.*-^'' — ReS. In leading order in A this becomes 

ImG ^ ~X"\lo - h^°'^)B'^"'\Lu - h^°'^) + 0{X''"'+^'^) . (A4) 

The negative semi-definitenes of the 1. h. s. o f (|A4| ) implies the positive semi-definiteness of B for values ofui which are 
no eigenvalues of h^°\ Assuming continuitytil for B^"^\u!) the positive semi-definiteness extends to all frequencies. 
Addressing the sign of the imaginary part of the full Green function we state tha t B > implies that there is a 



matrix V B which is real and symmetric as well. Defining D := V B A v B allows (A3) to be written as 



lmG = -A-^VB{l + D^) ^VbA-\ (A5) 

The expression in parentheses is manifestly positive semi-definite, thus its inverse as well. Since A~^^/B and ^/B A~^ 
are transposed to each other, the r. h. s. in (|A5D is negative semi-definite as a whole. This concludes the argument. 
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FIG. 19. Diagrammatic basis element of the terms in the sum in fig. |l^ 



Note that the momentum of the upper propagator cancels the momentum in the lower propagator up to multiples 
of Q = (tTjTTjTt)^, since we are interested in the average conductivity. At the end-vertices no momentum is added or 
subtracted. The case k = k' = q = q' can be combined with the case k = q = k' + Q = q' + Q. In both cases the 
momentum through the whole segment is zero. Using the elements from (0) one obtains 



2UT 



E 



UJ I, — ijJ X—ijJ n 



Nc,o{e)de , 



(Bl) 



where the short-hand w^/\ :~ ii-u^/x — T,{iLU^/x) is employed again. A similar expression is obtained in the combined 
cases k = k' = q+ Q = q' + Q and k = q + Q = k' + Q = q' for which the momentum through the segment is Q 



A. 



2UT 



E 



(B2) 



The segments which change the momentum are very important. The cases k = k' = q + Q = q' and k = q + Q = 
k' -f- Q = q' -I- Q yield together 



A, = - 



2UT 
Z 



E 



LOii—UJX—UJr, 



WiyA{iuj\) - w\A{iuJiy) 



-Nc,o{e)de . 



[wl - (7£)2 - A^iiu,)) (wl - (7e)2 - A^iicux)) 
The possible combinations of these cases are naturally generated by powers of the matrix 

^3 A2 

The sum of all these powers is a geometric series yielding finally 



(B3) 



2 
U 



(B4) 



(B5) 



Taking the (1,1) element ensures that the average current-current correlation function is calculated. The prefactor 
compensates the fact that the end vertices do not have the factor which is incorporated in A. Eq. ( pl| ) is a direct 

consequence of (B5). The representation ( ^2|) follow from (Bl), (B2), and (B3) by partial fraction decomposition and 
integration over the conductivity spectral density using the conductivity Green functions. 

Now we turn to (^). For x''''(0) we need only Ai at ium = since A^ vanishes at ium — and the factor 1 — A2 
then drops out. For Ai one obtains 
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2UT 
Z 

2UT 
Z 

2UT 



E 
E 

z S 



{w^+-fe)^ - A^(w^) 

K_(^£)2_A2(ZC.,.)) 

7V,,o(£)- 5(eiV,,o(e))/5£ 



2Afc,o(£)(^£ 



-(7e)2-A2(iw,) 



de 



-No{e)de 



00 



(w-E(w))2- (7e)2- A2(w) 



(B6a) 
(B6b) 
(B6c) 
(B6d) 
(B6e) 



The step to (B6b ) follow s from the fact that Ncfl{s) is even. Partial integration of the last term of the numerator 
leads to (B6c). Eq. (B6d) is obtained using ( p9| ) which emphasizes the importance of this relation between DOS and 
conductivity spectral density. In the end one obtains 



7 
U{T) 

= -(1/7-1). 

One obtains by substituting ( B7c ) in the denominator of (|l|) and comparing the resulting expression with ( |B7lj ). 
This completes the analytic derivation of the f-sum rule in the self-consistent 1/d approximation which is considered 
here. 



(B7a) 

(B7b) 
(B7c) 
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